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ABSTRACT 

This theoretical study of the trajectories of charged particles 
supplements an experimental project of molecular-ionic rearrangements 
using the magnetic focusing properties of a beta-ray spectrometer. The 
experimentally measured magnetic field was analytically represented 
by a twelfth-order polynomial. This field is axially symmetric, but 
non-homogeneous otherwise. The particular particle of concern was 
the hydrogen ion, H 2 + . The trajectory of this particle was computed 
from a second-order differential equation, assuming values for the 
kinetic energy and initial angle of the particle, and the magnetizing 
current. The solution was obtained by numerical integration using a 
CDC 1604 digital computer. The distinctive feature of these calcu- 
lations, in contrast to those used normally for a beta-ray spectrometer, 
is the large scattering (initial) angle, about 45°. 



2 



TABLE OF CONTENTS 



Section 




Page 


1 . 


Introduction 


9 


2. 


Discussion 


9 


3. 


The Magnet and its Field 


10 


4. 


Theory of the Trajectory 


14 


5. 


Expression for the Vector Potential 


16 


6. 


The Computer Solution 


17 


7. 


Results 


21 


8. 


Acknowledgements 


22 


9. 


Bibliography 


24 


Appendices 






1 . 


Data on the Magnetic Field 


25 


2. 


Useful Formulae in Cylindrical Coordinates 


28 


3. 


Power Series Expansion to Obtain the Vector 
Potential 


29 


4. 


The Differential Equation Development 


33 


5. 


Development of the Vector Potential 


37 


6. 


Pertinent Machine Programs 


39 


7. 


Graphs of Various Trajectories 


66 



3 



LIST OF ILLUSTRATIONS 



Figure 




Page 


1 . 


The Magnet 


11 


2. 


The Magnetic Field 


13 


3. 


The Standard Trajectory 


67 


4. 


Trajectories Showing Variations in Energy 


68 


5. 


Trajectories Showing Variations in Magnetic 
Field 


69 


6. 


Trajectories Showing Variations in Initial Angle 


70 



5 



TABLE OF SYMBOLS 



A The vector potential 

B The magnetic flux density of the field 

B(i) The coefficients computed for the polynomial which describes 

the magnetic field on-axis 

E The kinetic energy of the rebound particle , a function of the 

scattering angle 

E 0 The energy of the incident particle 

Fi, Gi Functions of z which are described later in detail 

H 0 (z) The axial component of the field on axis 

H 1 The i-th derivative of the function H 0 (z) 

IB The magnetization current 

k = p/q 

L The intercept distance measured on the z-axis 

0 The azimuthal angle in cylindrical coordinates 

p The momentum of the rebound particle (H 2 + ) 

q The magnitude of the electronic charge 

r The radial distance in cylindrical coordinates 

6 The scattering angle measured in the laboratory frame 

z The axial distance in cylindrical coordinates 



7 



1 . Introduction 

This presentation is part of an experimental project presently being 
conducted at the U. S. Naval Postgraduate School to verify a classical 
theory developed by Bates, Cook, and Smith [l ] of the capture of a 
light ion or atom from a target system by a fast projectile. 

Specifically, the theory involves a reaction of the type: 

H + + CH 4 - H + H + CH 3 

and is essentially a sequence of binary collisions. The process 
ensues with the bombardment of methane gas by low energy protons, 
whereby the resulting rebound proton and escaping hydrogen ion join to 
form the product H 2 + . The angle of departure (or scatter) is shown to 
be in the vicinity of 45 degrees, and is labeled 0. 

Since the expected cross section will be quite low (^1CT 18 cm 3 ), 
a detector could not be utilized, as is standard procedure, by placement 
at the scattering angle. [5 ] 

Therefore, the beta-ray spectrometer was chosen due to its prop- 
erty which will cause all H 2 + ions of the same departure angle and 
energy to converge. This study determines that point of convergence, 
or focus . 

2. Discussion 

It is known from elementary electromegnetic theory that a magnetic 
force, q{v x B], is produced on a charged particle moving at an angle 
to the direction of the field. This force causes the particle to sustain 
a change in velocity (direction only). The supposition is, for this 
field, that a particle which begins its motion on the axis of the field 
and at an angle with that field axis, will travel a curvilinear path and 
ultimately intersect the axis of the field for appropriate values of the 
field strength, energy and initial angle of the particle. This can be 
readily shown for a homogeneous field. 

The point at which the particle returns to the axis is called the 
intercept point, and the distance from the origin of the trajectory to 
the intercept point is the intercept distance (L) . This is where the 
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detector can be placed to analyze the particulate beam . The experimental 
usefulness is that if the intercept distance is known along with the 
field strength, and the energy of beam E 0 , then the initial angle may be 
evaluated. 

The theoretical problem, which is herein shown, is the evaluation 
of the intercept distance assuming values for the scattering angle. 

Thus, the following intermediate problems arose: 

A. A means to express the field, or more precisely, a continuous 
function which would represent the field for all points in the area of 
interest. 

B. A set of equations to specify the motion of the particle to 
incorporate the energy and charge of the particle, as well as the inter- 
action of the magnetic field. 

C. A numerical method to solve a differential equation which 
would result from the equations in (B) above. 

D. Interpretation and extrapolation from sets of trajectories to 
obtain these derivatives: 

— The change in intercept distance with 
respect to changes in energy 



-- The change in intercept distance with 
respect to changes in angle 



-- The change in intercept distance with 

respect to changes in magnetization current 

3. The Magnet and its Field 

The magnet which will be used for this project was previously in 
the USNPGS beta -ray spectrometer. It is essentially composed of a 
normal section of a right circular cylinder, and wound with cylindrically 
concentric layers of copper wire. The outer disks and cylindrical 
center section which form the container for the layers of wire are 
copper as are the cooling water tubes. The figure (Fig. 1) on the 
next page shows the magnet with its dimensions . 
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The experimental field strength of this magnet on-axis, H 0/ appears 
as a Gaussian curve as shown in Fig. 2. The quantities in Appendix 1 
for the values of the field are shown in the positive z direction. The 
maximum intensity of the field on-axis is at the geometrical center of 
the magnet, at the point formed by the intersection of lines ii and jj. 

B has cylindrical symmetry about the line jj , and is symmetrical about 
the center plane. 

In order to more easily discuss the field, the system of cylindrical 
coordinates should be established at this point. The origin is defined 
to be at the intersection of the lines ii and jj . The line jj will be 
the z axis, and the radial distance, r, will be measured perpendicularly 
from it. The angular component, 0, is as in (r, 0, z) , a right hand 
coordinate system. 

Measurements were made on this field for various values of field 
current by P. J. Kelly in May 1965 with an axial Hall Probe, and this 
data has been compiled in Appendix 1 . The axial component of the 
field intensity was measured initially on the primary axis of the magnet 
(z axfs) for a distance of 66 cm. on either side of the origin at 4 cm. 
intervals. This was done for various current settings between 2 and 
14 amperes . 

Also, the axial component of the field was measured off axis at 
four values of z, and for 8 values of 0, for each value of z. The 
radius r was varied from zero to seven cm. 

The conclusions concerning the field were that: 

(1) The field strength is linear with current over the entire 
range of current settings within an accuracy of ±0.1%. 

(2) The field is symmetrical on either side of the z = 0 
plane to ±0.01%. 

(3) The field is azimuthally symmetric to ±0.01%. 

To find the analytic expression for the field on axis, H 0 (z), where 
z assumes positive values only, a polynomial was fit to the data by 
the method of least-squares. This data is contained in Appendix 1 . A 
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twelfth order polynomial was chosen, and provided sufficient precision 
in the form: 



H 0 (z) = B l + B 2 z + B 3 z s + B 4 z 3 + B 6 z 4 + . . . + B 13 z 13 . (3-1) 



The thirteen coefficients were determined by the computer and a USNPGS 
subroutine LSQPOL, which is in PROGRAM CURVE included in Appendix 

6 . 

At this point, a comparison was made of the experimental field 
values and those values which would result from the polynomial field 
expression . 

To make this check of the values of the field off-axis, the following 
steps were taken: 

(1) The vector potential was computed for all values of r and z 
using the derived polynomial. 

(2) The axial component tof the field was derived, using the 
vector potential, and evaluated for various values of r and z. 

This computation is shown in detail in Appendix 1 . The agreement be- 
tween experimental and theoretical values was again quite good. 

4. Theory of the Trajectory 

The H 2 + ion has a positive charge eq^al in magnitude to that of the 
electron, and mass essentially equal to twice the mass of a single 
proton. The mass of the single electron is ignored. 

The path of this particle can be represented by this equation: 



[l + (r") s ] 



{k 3 - A 3 ] - r'(A) 




(4-1) 



where A is the vector potential, and k is defined as (mv/q). The 
prime above indicate differentiation with respect to z . A discussion 
of A, its derivatives, and k will follow. 

This equation for any charged particle can be derived from the 
equations of motion for that particle which begins its path on the 
primary axis of the magnetic field. This equation is derived in Appendix 
4. 
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The energy of the H 2 + particle after scattering is E, and can be 
shown to be: 



where E 0 is the kinetic energy of the particle which is incident on the 
target methane. The angle 8 is again the angle of scatter. 

A description of the collision process is as follows: 

The incoming proton strikes one of the hydrogens in the methane molecule. 
This is the first of the two binary collisions, and the only one which is 
assumed to be totally elastic. The initial particle moves off at an 
angle, called the departure angle; and the second particle (essentially 
a proton and an electron) moves towards the carbon atom at an angle 
equal to 90° - 8 . Since the carbon atom is assumed to be infinitely 
massive, compared to the hydrogen atom approaching it, this second 
particle is redirected in the second binary collision in the direction of 
the recoiling proton. 

A schematic of the collision process is as follows: 



E = 2 E 0 cos s 8 



(4-2) 




carbon 



In the laboratory frame, the initial energy of the system is: 



E 0 = 4 mj (v"i ) 2 



(4-3) 



and the final energy is: 



E = (2m) (v') 2 



(4-4) 



Conservation of momentum yields, for the first collision: 



P x = mv 



mv' cos 8 + mv" sin 8 



(4-5) 
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and, 



p y = mv 1 sin 6- mv' 1 cos 6=0 



(4-6) 



The last equation results in 



v“" = v 1 tan 6 



(4-7) 



and substituting into the equation (4^5) 



v" 1 = v* cos 6 . 



(4-8) 



With this relationship between the velocities , the final energy after 
the two collisions becomes 



If both collisions are elastic in nature, then it is quite easy to show 
that 6 must be 45° . The y-component of momentum produces the 
equation 



at the first collision. And it can only be satisfied when 

sin 6= cos 6 
or, when 6 = 45°. 

5. Expression for the Vector Potential 

The fact that the field is cylindrically symmetric carries with it 
the implication that there is no azimuthal component, since the curl 
B = zero. This field can be described by a vector potential which is 
always in the azimuthal direction. 

Using the calculated expression for magnetic field on-axis, H 0 (z), 
the vector potential expansion 



E = *(2m) (v‘) 2 



(4-4) 



= ^ (2m) v 3 cos 2 6 



= 2E 0 cos 2 6 
E = 2E 0 cos 2 6 



(4-2) 



Py = mv' sin 6- mv' cos 6=0 



(4-9) 




+ . . . + 




H 0 (z) + 
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is shown to be for any point (r,z): 



A = — [G t ] - — [G 2 ] +_L_[G 3 ] - — E — [G 4 ] + £ 

2 16 384 18432 1474560 



[g 6 ] 



r 11 



176947200 



[G 6 ] + 



29727129600 



[G 7 ] 



where 



H 0 (z) is the nth partial derivative of the function 



H 0 (z) with respect to z, and where G x = H 0 (z), and G 2 



w 



H 0 (z), etc. 



6. The Computer Solution 

It remains, at this point, to find the solution of the differential 
equation. The method chosen for this task was a process of numerical 
integration employing a digital computer. The subroutine INTEG-1 of the 
USNPGS Computer Library was selected because of its adaptability to 
this particular problem. The program in which INTEG-1 is contained 
is called PROGRAM ORBITS, and is listed as Program C in Appendix 6. 

The first four statements of this program (ORBITS) are standard 
according to the instructions for use of INTEG-1 . The three comment 
cards that follow describe ENER, OKAY, and P, which are names given 
various quantities used in the solution. 

B (1 ) through B (1 3) are the coefficients used in the polynomial 
as obtained by PROGRAM CURVE. The functions F(l) through F(6), 
and G(l) through G(7) are explained in Appendix 5. 

VECDZ is the symbol for the partial derivative of the vector 
potential with respect to z. VECDR represents the partial derivative 
of the vector potential with respect to r. VECPT is the vector 
potential . 

A second order differential equation can be reduced to a system 
of two first differential equations. Applying this technique to 
equation (4-1) yields 



_ fA dA i . [a bA 1 r a dA 

Cr' ]' = r' LA 77 J + (r') 3 LA 77 J - (r') 2 LA 7T 



[k 2 - A 3 ] 



a z 
[k 2 - A 3 ] 



] [A 



aA 

ar 



] 



[k 3 - A 3 ] 



[k 2 - A 2 ] 
(6-1) 
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[r] 



( 6 - 2 ) 



" = r - 

Letting 

x(l) - r and x(2) = r* , 

the equations (6-1) and (6-2) become 

[x(2)]' = [x(2) ] FUNC4 + [x(2)] 3 FUNC4 - [x(2)] 2 FUNC5 - FUNC5 
[x(l)]‘ = x(2) 
where 

FUNC4 = ^ 

o z 

[k 2 - A 2 ] 

FUNC5 = — 17 - 

or 

[k 2 - A 2 ] 

Further, FUNC4 and FUNC5 can be written in terms of FUNC1 , 
FUNC2 and FUNC3 as follows: 

FUNCl = k 2 - A 2 

FUNC2 = A AA 
dz 

FUNC3 = A AA 
dr 

and 

FUNC4 = 

FUNCl 

FUNC5 = FUNC3 
FUNCl 

OKAY is obtained in the following way. This is the machine 
language expression for k 2 where the energy can be entered in eV. In 
the MKS system, k must have the unit of webers/meter. This is the 
same as A, since B has the units of webers/meter 2 , and B = 7xA; 

k = m v in units of (kilogram) (meter/ sec) or webers/meter, 
q (coulombs) 
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It would be desirable to have k expressed in terms of the energy 
of the particle, E, and in the units of eV, as a great savings in con- 
version would be realized. 

Using the non-relativistic expression for the energy E, the 
following is obtained: 

E =_e 1_ p = (2mE) i 

2m 

But, since k 2 is needed 

k 3 - — 2? — = if E‘ is expressed in eV. 

q 2 q 2 q 



And, k 2 becomes 



[ 2 (ma s s of H 2 + ) 1 
charge of electron! 



E' 



Inserting the values of the constants, k 2 becomes the quantity OKAY: 

-f 



OKAY E f 4(1. 6725 x IQ' 87 \ 
1.6021 x 10 -19 j 



= {4.17576 x IQ - ® 3 E' 



or 

OKAY = (4.17576E - 08)*(ENER), as expressed in the program. 

The method of solution employed by INTEG-1 is a fourth order 
Runge-Kutta numerical integration approximation of the ordinary differ- 
ential equations. 

The computation is performed in the MKS system of units . 
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TO USE PROGRAM ORBITS (the main trajectory program) 

This program is written in the language of FORTRAN-60. Its 
successful employment depends on the correct use of the data cards, 
which are the last eight cards of the program itself. These eight cards 
contain the input information, and the quantities desired as output. 

An example of the data cards may be seen on the last page of Program 
C in Appendix 6. 

Data card #1 : This card is used for identification, and columns 2-32 
are available for this use. This title appears on all data output, as 
well as graphs . 

Data card #2 : This card contains the number of "runs", or trajectories, 
to be processed with a maximum of 9. This number is entered in 
column 1 of the card. Note, however, that only one graph is possible 
for the entire program . 

Data card #3 : The number of coupled first-order differential equations 
is specified by this card (maximum of 30). In this case, there are two 
and this number is posted, right justified, in columns 1 and 2, as 02. 
Data card #4 : The initial and final values of the independent variable, 
z, over which the integration process is to be carried out. The inte- 
gration step size is also stipulated. The integration can be processed 
in up to 3 segments, each with a different step size, thus: 

ZI -S- 1 » Zl - S2 «-Z2 S3 „ ZF„ 
or ZI— SS-»Z1- S 2 »-ZF , 
or ZI- S1 »ZF , 

where ZI refers to the initial value of z (in this case equal to zero), 

ZF refers to the final value of z, Zi refers to intermediate values of 
z, and Si are the step sizes. All of these numbers are entered in 
meters. The corresponding values, with decimal points, are placed 
in the above order in columns 1-10, 1 1-20, 21-30, etc. 

Data card #5 : Columns 1-10 contain the value of the energy E in 
electron-volts, again, with the decimal point. Columns 11-20 
contain the quantity P which adjusts the magnetic field according to 
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the base of 10 amperes magnetization current. For example, for a field 
of 7 amps, the quantity 0.70 would be entered anywhere in columns 
11 - 20 . 

Data card #6 : Here are the values of the initial conditions. The format 
is the same as for cards 4 and 5. The first value expressed is the 
initial value of r (equal to zero), and the next value is the slope of 

the function at that point, ( dr \ # This derivative may be considered 

\dz) 

to be also the tangent of the scattering angle, 8 . 

Data card #7 : This card controls data print-out. A blank card will 
suppress this output. Again, using the ten column sequence, this 
format is as follows. The first 8 columns contain the title for each 
variable output desired, and the last two columns contain the variable 
identifier. The variable identifiers are 00 for z, 01 for r, and 02 for 
r‘ . An example: Z DIST 00R DIST 01SLOPE 02 
Data card #8 : This last data card controls graph output. A blank 
card will suppress this output. Remember: that only ONE graph for 
each submitted program can be drawn and it corresponds to the first 
set of data. The information is contained in columns 1-20. Columns 
1-16 contain the graph title, and columns 17-20 indicate the variable 
identifiers. Columns 17-20 must always be 0100. The origin, 
direction of plot, and scale distances have already been set within 
the main program . 

Running Time : The time required for the computer solution is dependent 
on the step size and the total range of the independent variable. This 
running time has averaged about three minutes for a single trajectory. 

7. Results 

The outcome of this project was the creation of a single-package 
system which generates hydrogen ion trajectories in an axially 
symmetric field. 

Approximately 100 production trajectories were computed for 
values of energy E ranging from 97 eV to 450 eV; magnetizing current 
from 7 to 15 amperes; and scattering angle of 35-55 degrees. The 
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numerical output data and the graph plots showed these orbit con- 
figurations to be quite smoothly continuous and consistent. 

The error in evaluating the intercept distance L was due mainly 
in the uncertainty of reading the graph or interpolating the data print- 
out. These values of L ranged from 20 centimeters to about 80 centi- 
meters, with an estimated error of ±1.0%. The maximum radial distance 
of the orbits was around 8 cm. 

Comparison of these computer trajectories with those obtained by 
the fine wire technique used by Kelly showed only fair agreement to 
within ±10.0%. The disparities in the different sets of trajectories 
have not as yet been determined. 

This computer system was checked with a constant uniform 
(homogeneous) magnetic field against an analytic solution. The com- 
parative results were excellent (to within ±0.01%). 

Based on a Standard Trajectory of E 0 = 300 eV # 

I B = Magnet current = 12 amperes, 

0 = 45 degrees : 

= - 3 . 3 ± 0 . 2 cm /deg 
B 

-, T i = -9.7 ±0.8 cm/amp 

M B/ 0, E 0 

„ = 0.20 ±0.01 cm/eV 

& E oy Ig 

The above derivatives show that as 0 and Ig are increased, L is 
decreased. However, an increase in E 0 produces an increase in L. 

Each of the partial derivatives were obtained with respect to a change 
in a single parameter, with the other two parameters held constant. 

8. Acknowledgements 

In this limited space and manner, I note with pleasure all those 
people from whom I have received impetus towards the completion of 
this study. To them, I render my personal appreciation. 



bL 

be 



22 



I wish to thank Professors F. D. Faulkner and H. M. Martinez 
of the Department of Mathematics, Professors F. W. Terman and J. R. 
Ward of the Electrical Engineering Department, and Professors O. 

Heinz and R. L. Armstead of the Physics Department for their inspiration, 
encouragement and assistance. 

Prof. Ward provided invaluable assistance and advice as to the 
modification of the integration and draw systems of the main trajectory 
program. Prof. Heinz, at whose suggestion and under whose guidance 
the central project is being undertaken, contributed his ideas and 
enthusiasm. Prof. Armstead, as my thesis advisor, aided with his 
timely counsel, helpful suggestions and continued confidence; and to 
him, I am sincerely indebted. 

To all the personnel of the Computer Facility, USNPGS, and, in 
particular, Mrs. Robin Ekstrum, I salute for their extreme consider- 
ation. Without their cooperation, the task would have been manyfold 
more difficult. 

Finally, I am grateful to my wife, Dorothy, for her understanding 
and devotion throughout this undertaking, and for typing the draft copy. 



23 



BIBLIOGRAPHY 



1. Bates, D. R., Cook, C. J., and Smith, F. J. Classical Theory 
of Ion-Molecule Rearrangement Collisions at High Impact 
Energies, Proceedings of the Physical Society , v. 83, 1964: 49-57. 

2. Deutsch, M., Elliott, L. G. and Evans, R. D. Theory, Design, 
and Applications of a Short Magnetic Lens Electron Spectrometer. 
The Review of Scientific Instruments ., v. 15, no. 7, July 1944: 
178-195 

3. Hedgran, A. Precision Measurements of Nuclear Gamma-Radiation 
by Techniques of Beta -Spectros copy . Arkiv For Fysik , v. 5, no. 1, 
May 1951: 1-27. 

4. Keller, J. M., Koenigsberg, E. and Paskin, A. Ring Focus in 
a Thin Magnetic-Lens Beta-Ray Spectrometer. The Review of 
Scientific Instruments , v. 21, no. 8, August 1950: 713-721. 

5. Kelly, P. J. Analogue Measurement of Charged Particle Tra- 
jectories in an Inhomogeneous Magnetic Field. Thesis (MS) 
USNPGS, Monterey, California, December 1965. 

6. Rosell, F. E., Jr. Thin Lens Beta Spectrometer. Thesis (MS) 
USNPGS , Monterey, California , May 1956. 

7. Siegbahn, K. M. Beta- and Gamma-Ray Spectroscopy. Inter- 
science Publishers, 1955 . 

8. Verster, N. F. Theory of a Magnetic Lens Type Beta Ray Spec- 
trometer. Applied Scientific Research , v. B-l, 1950: 363-378. 



24 



APPENDIX 1 



DATA ON THE MAGNETIC FIELD 

The magnet which provides this inhomogeneous field was formerly 
employed in the beta-ray spectrometer which was reconstructed here 
at the USNPGS in the period 1955-5 6. Although much was known about 
the magnet itself, there existed no data on the magnetic field until 
measurements were made by Kelly [5 Jin 1965. 

This data has been compiled and is herein available for the field 
strength both on- and off-axis; and also, for the theoretically computed 
as well as experimentally measured values in Gauss. 

In the diagram of the Magnetic Field Strength in this appendix, 
there are two values listed for each (r, z) . The experimental values 
in the diagram are the upper ones listed in each square. 

The theoretically computed values are listed in squares and are 
the lower values . The method by which these values were computed 
is as follows: 

The field on-axis is given by: 

H 0 (z) = B l + B 3 z + B 3 z 2 + B 4 z 3 + B s z 4 + . . . + B l3 z 12 (3-1) 



(z expressed in meters) 



where: 








B( 


1) 


= 


8 .0498907E-02 


B( 


2) 


= 


2.3351 83 IE-02 


B( 


3) 


= 


-4 . 8664438E+00 


B ( 


4) 


= 


4 . 2887590E+01 


B( 


5) 


= 


-4 . 36261 83E+02 


B ( 


6) 


= 


3.9271101E+03 


B( 


7) 


= 


-2.2648265E+04 


B ( 


8) 


= 


8.2393030E+04 


B( 


9) 


= 


-1 .9359723E+05 


B (1 0) 


= 


2 . 9437706E+05 


B (1 1 ) 


= 


-2 . 8025472E+05 


B (1 2 ) 


= 


1 . 5200646E+05 


B (X 3) 


= 


-3.5865075E+04 



Since A ^ is the only non-zero component of the vector potential, the 
axial component of the field strength, H z , can be computed and compared 
to values obtained by the Hall Probe. 
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H z becomes 



Mf;( Cr *4 



where A^ is equation (5-1). 

And, after performing the indicated operations. 



H z (r,z) - G x - — — G 2 +-^— - G ; 



64 



2304 



G 4 + .... 



This is demonstrated in Appendix 6. 



( 11 ) 
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Upper values - experimental 
Lower values - theoretical 



APPENDIX 2 

USEFUL FORMULAE IN CYLINDRICAL COORDINATES 



Conversion from"] 


x=r cos 0 


Rectangular r 


y=r sin 0 


Coordinates J 


z=z 


The Gradient 




7 * = + I 


(Ul\ 0 + (**-) k 
\b<t>) \bz J 


Ur J r 


The Divergence 




S - A = 1 (S(rAr)'l 


m 1 (SA 0V f dAz 1 


r 1 Sr J 


r l a<S / l 3z J 


The Curl 




7x~A = 1 r.BAz, - 


■»A0)1 J 


r L $0 


dz J 


r dAz 


1 

<v 1°^ 

N 5? 

1 1 


L ar 




+ 1 TS (rA tf> ) aAr 




r L dr d 0 



The Laplacian 



7 2 $ = 
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APPENDIX 3 

POWER SERIES EXPANSION TO OBTAIN THE VECTOR POTENTIAL 
The well known expansion, equation (5-1) will now be derived. 
This expression describes the vector potential for an inhomogeneous 
magnetic field, which acts only in the azimuthal direction. The field 
on-axis, H 0 (z), and its even derivatives can be computed. 

Basically, the derivation of this expansion is as follows: 

Since the field is constant and involves no currents in the region of 
the field, 

7 x B = 0 . 

And, 

B = 7 x A, 

so that the vector equation for which a solution is desired is: 

7* x ( 7 x A) = 0 . 

Conditions for this field are: 

B 0 ~ 0, 

and A r = A z = 0. 

In cylindrical coordinates, 7 x B becomes 



7 X *B 


= I 


'6B Z d(rB q) 


* 


3B Z 


d B r " 


0 + 


1 


) 


dB r 




r 


_ d 0 d z _ 


T 


. ^ r 


dz _ 


r 


dr 


d 0 


each component must be equal to zero. 










dB z 


= d% 
















dr 


d z 














a) 


B r = 


i I 


SA Z 1 _ f d (rA 


<h )' 


= 1 


f _ d (rA 


0 ) 


1 


(2) 


r 1 


d 0 J L 




r 


1 


dz 




; 


B 0 = 


5A Z 


i 

cv 

> 

III 


0 














dr 


dz 
















B z = 


i r 


B (rA ^ )"]_ 3A r 


= 1 


f 3 < rA <p 


) 


4 


(3) 




r L 


dr J d 0 


r 


l 3r 




J 
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N> 



From equations 1 , 2 and 3 , 



£[H-Ml--h[j{- it (tA, j 



= 0 , and 



r[» * f] 



a 

dz 

dr 



dr A d A 

dz r dz 






and the partial derivative is zero as the coordinates are indepen- 

d z 



dent . 



Thus , 



5 s A I SA. A a 2 A _ 

dr 2 r dr r 2 dz 2 

CD 

Here, a solution is assumed in the form E a n (z) r n where 

n=0 n 

the coefficients a n (z) are differentiable. To determine the coefficients, 



-7“ = £ n a n (z) r n 1 

or ' 



and 



= £ n(n-l)a n (z) r 11 2 



and 



dA 

dz 



= Z a n ' (z)r n 



and 

&= S V Wr» 

Substituting in, and multiplying through by r 2 
Z (n) (n-l)a n (z)r n + Zn a n (z)r n - Za n (z)r n + Za n "(z)r n+2 =0 

Z { n(n-l) + n - l}a n (z)r n = - Z a n "(z)r n+2 
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For n=0 / the left hand side of the last equation reduces to -a 0 (z), and 
there exists no n=l term. Therefore, 

-a 0 (z) + n f 2 {n(n-l) + n-l]a n (z)r n + n J^ a n "(z)r n+2 = 0 

-a 0 (z) + £ { (n+2)(n+l)+n+l } a n+2 (z)r n+2 + £ a n "(z)r n+2 =0 

n=0 n=0 

-a 0 (z) + L { (n+3)(n+l)} a n+2 (z)r n+2 + L a n M (z)r n+2 = 0 
n=0 n=0 



-a 0 (z) + f { [ (n+3)(n+l)]^, +2 (z) +a n "(z)}r n+2 = 0 

n=0 

The following must then be true, as for all power series expansions: 
a 0 (z) = 0 / hence 
(n+3) (n+l)a n+2 (z) + a n "(z) = 0 

aijj+2 (z) = ~ a n" 

(n+3)(n+l) 

Whence, a recursion relationship yields the important needed para- 
meters. Obviously, there can be no terms for n even. 

For n = 1 a 3 (z) = ~ a i"t z ) 

8 

n = 3 a 5 (z) = ~ a a , M ( z ) , etc. 

24 

As A(r,z) = L a n (z) r n 
n=0 

= a^zjr + a 3 (z)r 3 + a 5 (z)r 5 + a 7 (z)r 7 + . . . 

Making the substitutions, 

A =a 1 (z)r-a 1 "(z)-j£ + (^) (t) V 'W - • . • 

And, to evaluate a x (z), solve for B z using 7 x A, 

Bz = = 2a a( z ) -| a i"( z )r 2 + ^ a," " (z)r 4 - . . . 
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Knowing that B is equal to H 0 (z) for r = 0, hence, 



H 0 (z) = B z (r=0) = 2a 1 (z) | /. a x (z) = £ H 0 (z) 

and the expression for the Vector Potential is obtained as 
A(r,z) = y H 0 (z) H 0 " (z) + ^ H 0 ""(z)- ... 
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APPENDIX 4 



THE DIFFERENTIAL EQUATION DEVELOPMENT 
The essentials of this operation involve simply applying Newton's 
second law in cylindrical coordinates to the magnetic force acting on 
the particle. Finally, time and the azimuthal variables can be elim- 
inated by means of the two motion constants, energy and the canonically 
conjugate angular momentum. Since the particle will begin its tra- 
jectory on the z-axis, the canonically conjugate angular momentum is, 
in fact, zero . 

Applying Newton's 2nd law, 

F mag = q(vxB) = m (a) 

where 




In cylindrical coordinates , the velocity vector is 

.A - A . A 

v = rr + r00 + zz, 
and, the acceleration vector is 

a* = (r - r 0 2 ) r + (2r 0 + r 0 ) 0 + (z) z. 

Knowing the vector potential , the magnetic field B and the 
quantity v x B can be found. 
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Hence: 




Equating each orthogonal cylindrical component, the following is 
obtained , 



m (r - r i> 2 ) = q( 5 -r— rA <t ) , 
o r 



m (2 r 0 + r *0) = q 



[- r f?( rA 0 ) - 7 ir( cA 0 )] - 



m(z) = q( 0 ~t— rA 0). 
o z 



These three equations can be reduced further to: 



(mr) = mr 0 2 + q 0 (rA 0) , 



dt 



■(mr 2 J) + qrA 0) = 0 , 



( 4 ) 

( 5 ) 

( 6 ) 



and 

-J"(mz) = qr 0 (A 0). 

From the middle equation above, it can be shown that the quantity 
conserved is the canonically conjugate angular momentum about the 
symmetry axis . Thus, 



= mr 2 0 + qrA 0 , 
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Replacing the independent variable of time ( t ) with the distance on 
the z-axis (z) yields 

2 + xj (£). + r * [ 






This operation can be performed as follows: 




Jto. * ^ .A • a 

v=rr+r00+zz , 

V 2 = V • V = (f ) 2 + (r j) 2 + (z) 2 # 

and 



r 



dr_ 

dz 




( 8 ) 
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and: 



m -553- ( 5 )’ * 3J ^-I £ )j - rn \jf~ " * Srj ( 9 ) 

”‘ 5 > [ir<*>] ■ i- {tt} [%■ - *] • «»> 



Combining the remaining equations 8,9, and 10, and then allowing 

that 



and 



Pc s 0 



k = 



the following equation results: 



[1 -Mr') 2 ] 





where the dots throughout this discussion indicate derivatives with 
respect to time, and the primes indicate derivatives with respect to 
the variable z. [4] 
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APPENDIX 5 

DEVELOPMENT OF THE VECTOR POTENTIAL 
To adequately describe the field at every point, and for use in 
the general trajectory previously derived, equation (4-1), the following 
quantities are needed: 

A = The Vector Potential , 
dA 



dz 



= The first partial derivative of A with respect to z. 



and 



dA 

dr 



= The first partial derivative of A with respect to r. 



The vector potential, found by the power series expansion in Appendix 
3, is 



A 2 16 ^ G ^ + 384 18432 ^ G ^ + * * *’ 

+ r 13 



29727129600 



-[g 7 ] , 



where, for the purposes of programming, the G functions are defined 
as follows: 

G: = H 0 (z) 



G 2 = 



G 5 = 



= * 2 H gfeL 

dz 2 

- 

dz 4 
a 6 H n 



dz b 



G = H n 
G? ' B z 12 
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The first partial of A with respect to r becomes, 



dA 

dr 



i[G 1 ]-|f-[G s ] + f^[G 3 ]- 



7r s 



18432 



[G 4 ] + 



9r 8 



1474560 



[G 5 ] 



llr 10 



176047200 



-[g 6 ] 



Ur 13 



29727129600 



*[G 7 ] 



The first partial of A with respect to z is a little more tedious, but 
can be shown to be, 

if- ifw * is-tw- w 



where the F functions are defined as follows (again for programming 
convenience) , 



Fi E 

Fp = 



_ a Gi 



az 

ac 2 

az 



P _ dG 3 

^ a — 



f 4 s 
f 6 

f 6 - 



az 

aG 4 

az 

aG 6 



a z 

sg 6 

az 



f 7 = o. 
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APPENDIX 6 

PERTINENT MACHINE PROGRAMS 
The pertinent computer programs used in the solution of this 
problem are listed in this appendix in the following manner: 



Program A 
Program B 
Program C 



Least-squares polynomial generation by the PROGRAM 
CURVE, and contains 7 pages. 

Theoretical off-axis field check labeled PROGRAM 
WATCHOUT, and contains 2 pages. 

The main trajectory program labeled PROGRAM ORBITS, 
and contains 15 pages. 
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Reverse Blank 



PROGRAM CURVE 

DIMENSION X ( 1 8 ) * F2(18)* W(18)* Y(18)* DELY (18)* B ( 18 ) * SB(18), T( 
118)* ST ( 1 8 ) * C ( 1 8 ) * CT ( 1 8 ) » A(30*30) 

READ 1 (X( I ) * 1 = 1 *18 ) 
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83 PPXPP=PPXPP+P( I )*P( I ) 00000560 

90 PM ( I ) =PT 00000570 

T(K)=PPXF/PPXPP 00000580 

PMXPM=PXP 00000590 



PXP=PPXPP 00000600 
A(K*1)=-ALPHA*A<K-1*1 )-BETA*A( K-2. 1 ) 00000610 
A(K*K-1 )=A( K-l,K-2 )-A(K-l*K-l )* ALPHA 00000620 
A(K.K)=1.0 00000630 
I F ( K— 3 ) 150*150*110 00000640 
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1 ERRORS/) 00001990 

188 FORMAT (23H0 LEGENDRE POLYNOM IALS/45H COEFFICIENTS OF Y=C1*L1+C2*L00002000 
12+ETC AND ERRORS/) 00002010 

600 FORMAT (41H1 CO EFFICIENTS OF Y=B 1+B2*X+ETC AND ERRORS/) 00002020 

601 FORMAT ( 3H T ( I 2 * 2H ) = 1PE 1 5* 7 , 6H ERRT = E 10 • 3 > 3H T ( OP I 2 * 2H ) = 1 PE 1 5 . 7 , 6H 00002030 



1 ERRT = E 1 0* 3 * 3H T ( OP 1 2 , 2H ) = 1 PEI 5 . 7 , 6H ERRT = E10.3> 00002040 

602 FORMAT { 3H C ( I 2 * 2H ) = 1PE 1 5. 7 * 6H E RRC = E10.3.3H C ( OP I 2 , 2H ) = 1 PE 1 5 . 7 > 6H 00002050 

1 ERRC = E 1 0 • 3 * 3H C ( OP 1 2 ♦ 2H ) = 1 PEI 5 . 7 * 6H ERRC=E10.3> 00002060 

603 FORMAT( I6*1P5E16#7) 00002070 

END 00002080 
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PROGRAM ORBITS 

DIMENSION X ( 30 ) *XDOT (30)»C(15)*B(15).F(15)»G(15) 
C<10)=1 

CALL INTEG1(T*X * XDOT » C ) 
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SUBROUTINE USES FOURTH-ORDER RUNGE-KUTTA METHOD TO INTEGRATE 00000020 



ORDINARY DIFFERENTIAL EQUATIONS SUPPLIED BY THE CALLING 00000030 
PROGRAM. SINCE THESE EQUATIONS MUST BE REPEATEDLY QUOTE OOOOOOAO 
CALLED UNQUOTE BY THIS PROGRAM, A COMPUTED GO TO IS HERE 00000050 
USED TO ENSURE RETURN TO THE CORRECT PART OF THIS SUBROUTINE. 00000060 
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15 < = K + 1 00000750 

16 CONTINUE 00000760 

IF ( J - 1)17,18,19 00000770 

17 PRINT 209 00000780 

GO TO 423 00000790 
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45 IF(XMODF(NOPTS, I NCPR ) ) 5 4 , 48 , 54 00002140 

46 IPAGE = IPAGE + 1 00002150 

I F ( NR - 1)1046,1047,1046 00002160 

1046 PRINT 218, (ITITLE(I), 1 = 1*6), IPAGE, ITITLE(7), ( JT I TLE ( I ) , I = 1 , 8 ) 00002 170 

PRINT 219 00002180 



GO TO 47 00002190 

1047 PRINT 1218* (ITITLE(I)* 1=1*6)* IPAGE* (JTITLE(I)* 1=1*8) 00002200 

PRINT 219 00002210 

47 PRINT 219 00002220 

218 FORMAT <1H1*///*20X*6A8*10X*5HPAGE *Il*14H OF OUTPUT FOR * A8 */////* 000022 30 
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70 CONTINUE 00002860 

GO TO 72 00002870 

71 PRINT 224 00002880 

224 FORMAT (//////*76H STOP WITH THE ABSOLUTE VALUE OF A DEPENDENT VAR00002890 

1 1 ABLE GREATER THAN 1.0E+12. */,57H INTEGRATION PROBABLY UNSTABLE. 00002900 



2 TRY A SMALLER STEP SIZE. ,26HNO GRAPHS WILL BE PLOTTED* ) 00002910 

GO TO 330 00002920 

72 DT = C ( 1 3 ) 00002930 

IF ( T I - TF) 73 *73*80 00002940 

73 I F ( T - TF ) 75 * 74 * 74 00002950 
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1977 SMLY = YMIN 00003540 

1970 CONTINUE 00003550 

TX(.l) = 0. 00003560 

TX ( 2 ) = 0. 00003570 

TX ( 3 ) = 0, 00003580 



oooooooooo oo 

o«-h 

l n'0'0'0'0'0'C' 0v0v0 h-h- 

co co co co co rn co co 

o o o o o oo o o o o o 
000000*0 000 oo 

oooooooooo oo 

o o O O O OTO o o o o o 

4 

f 
i 



o o O o o o 

m in vO x co 
h- x x x x x 
co co co ro co co 

o o o o o o 
o o o o o o 
O O 0,0 o o 
o o o o o o 



o o o 


o 


o o 


o 


o 


o 


o 


o 


CM CO <J- 


in 


x co 


O' 


o 


CM 


CO 


<t 


00 CO 00 


00 


00 00 


00 


O' 


O' 


O' 


O' 


CO CO CO 


CO 


CO CO 


CO 


CO 


CO 


CO 


co 


o o o 


o 


o o 


o 


o 


o 


o 


o 


o o o 


O 


o o 


o 


o 


o 


o 


o 


O O O 


o 


o o 


o 


o 


o 


o 


o 


o o o 


o 


o o 


o 


o 


o 


o 


o 



to 

< 



o 

o 

o 

CM 

o 

o 



x x x 
JOO 



to CD CO 
M II H 



X 

a 

Z • 

to o 

II II 



o o 

II II 



— X 
I— X X 
v: ^ * 
m 

II — CM 
II 



o 

o 



o 

e 

o 



o 


#» 




*> 








«— i 


o 


o 




o 


«. 


• 




• 


o 


* 




** 


e 


r-H 




« — i 


CM 


# 




• 


•V 


r. 




•* 


o 


a 




a 




a 




a 


o 


x 




X 










r-H 


X 




X 


o 


t— • 




*— 1 


• 


*. 




o 




a 




a 




a 




a 


• 


CD 




CD 


* 


< 




< 


UJ 


a 




a 


_J 


e> 




o 


X 


o 




o 




* 






X 


> 




> 




a 




a 




3 




3 


_l 


U 




u 


a 


c 




Q 


CD 


o 




O 


< 


z 




z 




* 








H 




CM 


o 


X 


a 


X 




•> 


a 


*" 


J »— 1 


1 — 1 


•> 


CM 




X 


CO 


X 


J X 




— a 





■ 4 * in X 

x x x 



CM CO 

X X 



<t in 
X X 



X — X 

5 II 

o> <— < < < 

ii ~ a > ii a 

UJLU OK Cl 

-J _J _J 3 a 
U I- K J U w J 
CD x •— i J C CO J 
< I- X < O < < 

JhmUI JU 



tO X •* 

I— l/) 00 
0 _ < O' 

Z — I 
D * *M 
I z o 

~ I CO 



CM 



CO 

to I- o 
i- to co 
a < co 
Z a 
ID * CO 
CM Z O 



UJ 

a 



_J 

UJ 

GO 

< 



> 

a 

X 

<3 

c 

o 

5 ; 



o 

o 



o 

o 



UJ 

CO 

< 



> 

a 

D 

u 

c 

o 

z 



o 

o 



a ii 

X > 

_) a 
D 3 
51 U 

— c 

a o 

— 51 

00 



X 

1 C II 4 3 
CM < 

co > ii a 
a c 
O X) -J 
1 — U IlI — I 
Q CD J 

o o < < 

O Z a U 

a m 

a cm 



d : 
z « 



■ > CO > 

i a a 



CO 


CO 


o 


X 


•> 


X X 




a 


* * ro 


CO 


CM 


<t X CO 


X 


CO 


X CO CO 




o 


o — -co o 


to 


x o 


to x ♦CM 


X 


to CO 


x to O CO 


a 


< CO 


a < o co 


z 


a 


z a o * 


3 


•> <r 


D * H CO 


CO Z 


o 


Z O — -CO 




1 CO 


•w a co 


X 




I z ~ 


<f 5 : 


a ii 


<f IS X 


< 


•— i 


< 1 


n a 


x > 


ii a 1 


Q 


a a 


o u 


a 


3 3 


a a a 


a a 


Z U 


a a z z 


co a 


— o 


cd a * 


< < 


a o 


< < a a 


a u 


— z 


a u — — 



X oo 
CM CM 
CO CO 



a 

CM 



o x 
CO CO 
CO CO 



62 



SS£S^°S 0000000000 °0000000 

^^Ci 0000 oo o°oohhhh H hhhhhS 

OOOOonnSn ^ OOOOOCOOOO °°000 

OOnonoSnSS 2 0000000000 °0000 

SSS22 OOOC 9 0 °00000000000000 

ooooooooo$oooooooooooo6oo§ 



o o 

r-H (\1 
<\J f\J 
<t <J- 

o o 
o o 
o o 
o o 



V o 

V tu 



o 

vD 



vO 

<\J 



cl 

z 



h- 

fVJ 



Z CL Z ( 
— O — ( 
cl t- cl i 
CL 10 CL t 



o a 




UJ 51 


o 


l- o 


CO 


UJ u 


* 


_J 




a z 




x uj 




O UJ 


< 


vj no 




Z UJ 


X •> 

o — 


UJ > 




UJ < 




CO X 


U H» 


to cl 
< o 


X U 


X u_ 


u — 


CL O 


1 - o 

CO 


o UJ 


•> w 


U 1 


U X 


_J 


a 


Q < 




UJ u 


h- • 


_J 


o - 


-J to 


o 


< z 


« ro 


u 3 


X — 


CL 


u 


z 


* X 


3 X 


1 — 


cl r- 


» 


co 




ti i #> 


z m 


z -« 


z -» 


O ~ 


— — o 


* 


(M U • 


UJ 


< 


X UJ 


1— * 


h- X 


h- ~ | 


I— 


3 o 


X 


iZ CO — 


ro X 


CL w o 


<f in 


X *-• 


•> * 


UJ «~ 


X X 


z z u 


X X 


t-H O 


•w w 


b- ~ u 


h~ h- 


3 lO 


< < 


ozu 




CL UJ 


CL CL O CO X 0 


O O Z D •- z 


u. U. UJ U) O K, 



o 

★ 



< 

* 



u 

+ 



-) 

'd- 

< 

+ 



< 

* 



V 

< 



X 

Q 

★ 



^ o m 

• • 

* o o 

- II II 

) — — 

■ t-< r\i 

* I- h- 

» u u 



in o 

• • 

O H 

II II 

CO <f 

I— h- 

u u 



Nj- I- 

► u 



z “> 

Z w 
X 



II + rH 



o u 

O K 



II II II 



o u 

O X 



CsJ 

CO 

CO 



CO 

CO 



o Z II Z o 

• Z Z — ' • 

VO h- 

•H -) r-i ~ 

»l II X II 

z “> • “> II z 

~ CL ~ ~ CL 

03<t*-«in~«o3 

") H hDQ 

w iij o y o ■ — ^ uj z z 

UCLQCQxUCLLULjJ 



CO <f 



in 



63 



G 

O 

(0 

G 

(0 

a 

x 

<D 



CD 

JG 



<d 

r a 



c 

c 



fC 5 

Q 

4 -> 

G 

a 

c 



to 

TD 

CD 

,G 

H 



d' 

O 



CL 

c 



a 

>1 

H 



O 

£ 



O 

a: 



o 

<r 



< 



iT\COD 
C- c\i O 1 — 1 
O • O Cl 

• i— i • <x 
O r-H £V 

o 

o c. i 

• o 

o o o 
c: c 



» C°l 



Cl 

NJ 



64 



Typical Output Data 
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STOP 



APPENDIX 7 

GRAPHS OF VARIOUS TRAJECTORIES 
This appendix contains these four figures which show trajectories 
as produced by PROGRAM ORBITS. 



Figure 3. 


The Standard Trajectory 


Figure 4. 


Variation in Energy 


Figure 5 . 


Variation in Magnetic Field 


Figure 6. 


Variation in Initial (scattering) Angle 
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6 7 



NOTE: From the origin in the lower right hand corner, the radial distance (r) is plotted versus 
the axial distance (z) , and both are in centimeters. Further, z is positive to the left. 
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The above curves range from 290 eV (the lowest curve) to 310 eV (the highest). 



o 




g § 
1 ” 



g c 

N <0 

n > 



00 




+1 



CO 

3 3 



■a 

•0 

> 

a> 



+1 



10 


C 

(0 






r 

o 

ts 


TJ 

Q> 


«3 

N 

13 

Q> 


S' § 
1 1 


0) 


Cm 


C 


® 10 


"(0 

u, 

H 


a> 

x: 

H 


o» 

(0 

£ 


5 2 

E - 6 



o> 

£ 



6g 



ie from 290 eV (the lowest curve) to 310 eV (the highest). 




69 




5 « 

1 I 



INITIAL DISTRIBUTION LIST 



No 



1 . Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 

2 . Library 

U. S. Naval Postgraduate School 
Monterey, California 93950 

3 . Defense Atomic Support Agency 
Department of Defense 
Washington, D. C. 

4. Prof. Robert L. Armstead 
Physics Department 
U. S. Naval Postgraduate School 
Monterey, California 93940 

5 . Prof . Otto Heinz 

Physics Department 

U. S. Naval Postgraduate School 

Monterey, California 93940 

6. Dr. Charles J. Cook 

Stanford Research Institute (A220) 
Menlo Park, California 

7. LtCol Edward Saunders 
Physics Department, USMA 
West Point, New York 

8. Capt. Patrick J. Kelly 
Canoe Hill Road 

New Canaan, Connecticut 

9. Lt. Thomas O. Bush 
Ordnance Engr. Curricula 

U. S. Naval Postgraduate School 
Monterey, California 93940 

10. Prof. Sam E. Gagliano 

Jefferson County Board of Education 
Birmingham, Alabama 

1 1 . Maj (Ret) F . J . Petry 
66 Walnut Lane 
Hicksville, New York 



71 



. Copie 
20 

2 

1 

1 

1 

1 

1 

1 

1 



12. Lt. Carl S. Gagliano 
Air Defense Board 

Fort Bliss , Texas 1 

13. Mr. Sam E. Gagliano , Jr. 

Boeing Co., Redstone Arsenal 

Huntsville, Alabama 1 

14. Capt. Ross A. Gagliano 
Student Detachment, OMR 

Fort Monmouth, New Jersey 2 



72 



Unclassified 

Security Classification 



DOCUMENT CONTROL DATA * R&D 

(Security claeelficatlon ol title, body ot abetted and Indexing annotation muet be entered when the overall report ie c lee allied) 



1 ORIGINATIN G ACTIVITY (Corporate author) 

U. S. Naval Postgraduate School 
Monterey, California 93940 



2a REPORT SECURITY CLASSIFICATION 

Unclassified 



Zb CROUP 



3. REPORT TITLE 

Theoretical Trajectories of Charged Particles in an Inhomogeneous Magnetic 
Field 



4 DESCRIPTIVE NOTES (Type ot report and Inctueive datee) 

Masters Thesis in Physics, May 1966 



5 A UTHO R(S) (Laet name, first name, initial) 



Gagliano , Ross A . 



Captain 



U. S. Army 




This theoretical study of the trajectories of charged particles supple- 
ments an experimental project of molecular-ionic rearrangements using the 
magnetic focusing properties of a beta-ray spectrometer. The experimentally 
measured magnetic field was analytically represented by a twelfth-order . 
polynomial. This field is axially symmetric, but non-homogeneous other- 
wise. The particular particle of concern was the hydrogen ion, H s + . The 
trajectory of this particle was computed from a second-order differential 
equation, assuming values for the kinetic energy and initial angle of the 
particle, and the magnetizing current. The solution was obtained by 
numerical integration using a CDC 1604 digital computer. The distinctive 
feature of these calculations, in contrast to those used normally for a 
beta-ray spectrometer, is the large scattering (initial) angle, about 45°. 



DD 



FORM 

1 JAN 04 



1473 



I 

73 Unclassified 

Security Classification 



Unclassified 

Security Classification 



14 . 


LINK A 


LINK B 


LINK C 




KOLB 


WT 


ROLE 


WT 


ROLE 


WT 


Molecular-Ionic Rearrangements 
Beta-Ray Spectrometer 
Hydrogen Ion Trajectories 
Numerical Integration 















INSTRUCTIONS 



1. ORIGINATING ACTIVITY: Enter the name and address 
of the contractor, aubcontractor, grantee. Department of De- 
fenae activity or other organization (corporate author ) issuing 
tha report. 

2a. REPORT SECUWTY CLASSIFICATION: Enter the over- 
all security classification of the report. Indicate whether 
"Restricted Data" ia included. Marking is to be in accord- 
ance with appropriate security regulations. 

2b. GROUP: Automatic downgrading ia specified in DoD Di- 
rective 5200.10 and Armed Forces Industrial Manual. Enter 
the group number. Also, when applicable, show that optional 
markings have been used for Group 3 and Group 4 as author- 
ized. 

3. REPORT TITLE: Enter the complete report title in all 
capital letters. Titles in all cases should be unclassified. 

If a meaningful title cannot be selected without classifica- 
tion, ahow title classification in all capitals in parenthesis 
immediately following the title. 

4. DESCRIPTIVE NOTES: If appropriate, enter the type of 
report, e.g., interim, progress, summary, annual, or final. 

Give the inclusive dates when a specific reporting period ia 
covered. 

5. AUTHOR(S): Enter the name<s) of author(a) *» shown on 
or in the report. Enter last name, first name, middle Initial. 

If military, show rank and branch of service. The name of 
the principal author is an abaolute minimum requirement. 

6. REPORT DATE: Enter the date of the report as day, 
month, year, or month, year. If more than one date appears 
on the report, use date of publication. 

7a. TOTAL NUMBER OF PAGES: The total page count 
should follow normal pagination procedures, Le., enter the 
number of pages containing information. 

76. NUMBER OF REFERENCES: Enter the total number of 
references cited in the report. 

8a. CONTRACT OR GRANT NUMBER: If appropriate, enter 
the applicable number of the contract or grant under which 
the report waa written. 

86, 8c, & 8 d. PROJECT NUMBER: Enter the appropriate 
military department identification, such as project number, 
subproject number, system numbers, task number, etc. 

9a. ORIGINATOR'S REPORT NUMBER(S): Enter the offi- 
cial report number by which the document will be identified 
and controlled by the originating activity. This number must 
be unique to this report. 

9b. OTHER REPORT NUMBER(S): If the report has been 
assigned any other report numbers (either by the originator 
or by tha aponaor), also enter this numbers). 

10. AVAILABILITY/LIMITATION NOTICES: Enter any lim- 
itations on further dissemination of the report, other than those 



DD i j?5“. 1473 (BACK) 



imposed by security classification, using standard statements 
such as: 

(1) "Qualified requesters may obtain copies of this 
report from DDC" 

(2) "Foreign announcement and dissemination of this 
report by DDC ia not authorized." 

(3) "U. S. Government agencies may obtain copies of 
this report directly from DDC. Other qualified DDC 
users shall request through 



(4) "U. S. military agencies may obtain copies of this 

report directly from DDC Other qualified users 
shall request through 



(5) "All distribution of this report is controlled. Qual- 
ified DDC users shall request through 



If the report has been furnished to the Office of Technical 
Services, Department of Commerce, for sale to the public, indi- 
cate this fact and enter the price, if known 

1L SUPPLEMENTARY NOTES: Use for additional explana- 
tory notes. 

12. SPONSORING MILITARY ACTIVITY: EiXer the name of 
the departmental project office or laboratory sponsoring (pay- 
ing for) the research and development Include address. 

13. ABSTRACT: Enter an abstract giving a brief and factual 
summary of the document indicative of the report, even though 
it may also appear elsewhere in the body of the technical re- 
port. If additional space is required, a continuation sheet ahall 
be attached. 

It is highly desirable that the abstract of classified reports 
be unclassified. Each paragraph of the abstract shall end with 
an indication of the military security classification of the in- 
formation in the paragraph, represented as (TS). (S). (C), or (U). 

There ia no limitation on the length of the abstract. How- 
ever, the suggested length is from 150 to 225 words. 

14- KEY WORDS: Key words are technically meaningful terms 
or short phrases that characterize a report and may be used as 
index entries for cataloging the report. Key words must be 
selected so that no security classification is required. Identi- 
fiers, such as equipment model designation, trade name, military 
project code name, geographic location, may be uaed as key 
words but will be followed by an indication of technical con- 
text. The assignment of links, roles, and weights ia optional. 



Unclassified 

Security Classification 



of charged part 







